Breakpoint characterization of a rare alpha0 -thalassemia deletion using targeted locus amplification on genomic DNA.

INTRODUCTION: The high-sequence homology of the α-globin-gene cluster is responsible for microhomology-mediated recombination events during meiosis, resulting in a high density of deletion breakpoints within a 10 kb region. Commonly used deletion detection methods, such as multiplex ligation-dependent probe amplification (MLPA) and Southern blot, cannot exactly define the breakpoints. This typically requires long-range PCR, which is not always successful. Targeted locus amplification (TLA) is a targeted enrichment method that can be used to sequence up to 70 kb of neighboring DNA sequences without prior knowledge about the target site. METHODS: Genomic DNA (gDNA) TLA is a technique that folds isolated DNA, ensuring that adjacent loci are in a close spatial proximity. Subsequent digestion and religation form DNA circles that are amplified using fragment-specific inverse primers, creating a library that is suitable for Illumina sequencing. RESULTS: Here, we describe the characterization of a rare 16 771 bp deletion, utilizing gDNA TLA with a single inverse PCR primer set on one end of the breakpoint. Primers for breakpoint PCR were designed to confirm the deletion breakpoints and were consequently used to characterize the same deletion in 10 additional carriers sharing comparable hematologic data and similar MLPA results. CONCLUSIONS: The gDNA TLA technology was successfully used to identify deletion breakpoints within the alpha-globin cluster. The deletion was described only once in an earlier study as the --gb , but as it was not registered correctly in the available databases, it was not initially recognized as such.

Human Bone Marrow-Resident Natural Killer Cells Have a Unique Transcriptional Profile and Resemble Resident Memory CD8+ T Cells.

Human lymphoid tissues harbor, in addition to CD56bright and CD56dim natural killer (NK) cells, a third NK cell population: CD69+CXCR6+ lymphoid tissue (lt)NK cells. The function and development of ltNK cells remain poorly understood. In this study, we performed RNA sequencing on the three NK cell populations derived from bone marrow (BM) and blood. In ltNK cells, 1,353 genes were differentially expressed compared to circulating NK cells. Several molecules involved in migration were downregulated in ltNK cells: S1PR1, SELPLG and CD62L. By flow cytometry we confirmed that the expression profile of adhesion molecules (CD49e-, CD29low, CD81high, CD62L-, CD11c-) and transcription factors (Eomeshigh, Tbetlow) of ltNK cells differed from their circulating counterparts. LtNK cells were characterized by enhanced expression of inhibitory receptors TIGIT and CD96 and low expression of DNAM1 and cytolytic molecules (GZMB, GZMH, GNLY). Their proliferative capacity was reduced compared to the circulating NK cells. By performing gene set enrichment analysis, we identified DUSP6 and EGR2 as potential regulators of the ltNK cell transcriptome. Remarkably, comparison of the ltNK cell transcriptome to the published human spleen-resident memory CD8+ T (Trm) cell transcriptome revealed an overlapping gene signature. Moreover, the phenotypic profile of ltNK cells resembled that of CD8+ Trm cells in BM. Together, we provide transcriptional and phenotypic data that clearly distinguish ltNK cells from both the CD56bright and CD56dim NK cells and substantiate the view that ltNK cells are tissue-resident cells, which are functionally restrained in killing and have low proliferative activity.

Comprehensive multi-center assessment of small RNA-seq methods for quantitative miRNA profiling.

RNA-seq is increasingly used for quantitative profiling of small RNAs (for example, microRNAs, piRNAs and snoRNAs) in diverse sample types, including isolated cells, tissues and cell-free biofluids. The accuracy and reproducibility of the currently used small RNA-seq library preparation methods have not been systematically tested. Here we report results obtained by a consortium of nine labs that independently sequenced reference, 'ground truth' samples of synthetic small RNAs and human plasma-derived RNA. We assessed three commercially available library preparation methods that use adapters of defined sequence and six methods using adapters with degenerate bases. Both protocol- and sequence-specific biases were identified, including biases that reduced the ability of small RNA-seq to accurately measure adenosine-to-inosine editing in microRNAs. We found that these biases were mitigated by library preparation methods that incorporate adapters with degenerate bases. MicroRNA relative quantification between samples using small RNA-seq was accurate and reproducible across laboratories and methods.

Brain Transcriptomic Analysis of Hereditary Cerebral Hemorrhage With Amyloidosis-Dutch Type.

Hereditary cerebral hemorrhage with amyloidosis-Dutch type (HCHWA-D) is an early onset hereditary form of cerebral amyloid angiopathy (CAA) caused by a point mutation resulting in an amino acid change (NP_000475.1:p.Glu693Gln) in the amyloid precursor protein (APP). Post-mortem frontal and occipital cortical brain tissue from nine patients and nine age-related controls was used for RNA sequencing to identify biological pathways affected in HCHWA-D. Although previous studies indicated that pathology is more severe in the occipital lobe in HCHWA-D compared to the frontal lobe, the current study showed similar changes in gene expression in frontal and occipital cortex and the two brain regions were pooled for further analysis. Significantly altered pathways were analyzed using gene set enrichment analysis (GSEA) on 2036 significantly differentially expressed genes. Main pathways over-represented by down-regulated genes were related to cellular aerobic respiration (including ATP synthesis and carbon metabolism) indicating a mitochondrial dysfunction. Principal up-regulated pathways were extracellular matrix (ECM)-receptor interaction and ECM proteoglycans in relation with an increase in the transforming growth factor beta (TGFß) signaling pathway. Comparison with the publicly available dataset from pre-symptomatic APP-E693Q transgenic mice identified overlap for the ECM-receptor interaction pathway, indicating that ECM modification is an early disease specific pathomechanism.

Immune stimuli shape the small non-coding transcriptome of extracellular vesicles released by dendritic cells.

The release and uptake of nano-sized extracellular vesicles (EV) is a highly conserved means of intercellular communication. The molecular composition of EV, and thereby their signaling function to target cells, is regulated by cellular activation and differentiation stimuli. EV are regarded as snapshots of cells and are, therefore, in the limelight as biomarkers for disease. Although research on EV-associated RNA has predominantly focused on microRNAs, the transcriptome of EV consists of multiple classes of small non-coding RNAs with potential gene-regulatory functions. It is not known whether environmental cues imposed on cells induce specific changes in a broad range of EV-associated RNA classes. Here, we investigated whether immune-activating or -suppressing stimuli imposed on primary dendritic cells affected the release of various small non-coding RNAs via EV. The small RNA transcriptomes of highly pure EV populations free from ribonucleoprotein particles were analyzed by RNA sequencing and RT-qPCR. Immune stimulus-specific changes were found in the miRNA, snoRNA, and Y-RNA content of EV from dendritic cells, whereas tRNA and snRNA levels were much less affected. Only part of the changes in EV-RNA content reflected changes in cellular RNA, which urges caution in interpreting EV as snapshots of cells. By comprehensive analysis of RNA obtained from highly purified EV, we demonstrate that multiple RNA classes contribute to genetic messages conveyed via EV. The identification of multiple RNA classes that display cell stimulation-dependent association with EV is the prelude to unraveling the function and biomarker potential of these EV-RNAs.

Detecting PKD1 variants in polycystic kidney disease patients by single-molecule long-read sequencing.

A genetic diagnosis of autosomal-dominant polycystic kidney disease (ADPKD) is challenging due to allelic heterogeneity, high GC content, and homology of the PKD1 gene with six pseudogenes. Short-read next-generation sequencing approaches, such as whole-genome sequencing and whole-exome sequencing, often fail at reliably characterizing complex regions such as PKD1. However, long-read single-molecule sequencing has been shown to be an alternative strategy that could overcome PKD1 complexities and discriminate between homologous regions of PKD1 and its pseudogenes. In this study, we present the increased power of resolution for complex regions using long-read sequencing to characterize a cohort of 19 patients with ADPKD. Our approach provided high sensitivity in identifying PKD1 pathogenic variants, diagnosing 94.7% of the patients. We show that reliable screening of ADPKD patients in a single test without interference of PKD1 homologous sequences, commonly introduced by residual amplification of PKD1 pseudogenes, by direct long-read sequencing is now possible. This strategy can be implemented in diagnostics and is highly suitable to sequence and resolve complex genomic regions that are of clinical relevance.

Flexible and Scalable Full-Length CYP2D6 Long Amplicon PacBio Sequencing.

Cytochrome P450 2D6 (CYP2D6) is among the most important genes involved in drug metabolism. Specific variants are associated with changes in the enzyme's amount and activity. Multiple technologies exist to determine these variants, like the AmpliChip CYP450 test, Taqman qPCR, or Second-Generation Sequencing, however, sequence homology between cytochrome P450 genes and pseudogene CYP2D7 impairs reliable CYP2D6 genotyping, and variant phasing cannot accurately be determined using these assays. To circumvent this, we sequenced CYP2D6 using the Pacific Biosciences RSII and obtained high-quality, full-length, phased CYP2D6 sequences, enabling accurate variant calling and haplotyping of the entire gene-locus including exonic, intronic, and upstream and downstream regions. Unphased diplotypes (Roche AmpliChip CYP450 test) were confirmed for 24 of the 25 samples, including gene duplications. Cases with gene deletions required additional specific assays to resolve. In total, 61 unique variants were detected, including variants that had not previously been associated with specific haplotypes. To further aid genomic analysis using standard reference sequences, we have established an LOVD-powered CYP2D6 gene-variant database, and added all reference haplotypes and data reported here. We conclude that our CYP2D6 genotyping approach produces reliable CYP2D6 diplotypes and reveals information about additional variants, including phasing and copy-number variation.

Full-Length Mitochondrial-DNA Sequencing on the PacBio RSII.

Conventional mitochondrial-DNA (MT DNA) sequencing approaches use Sanger sequencing of 20-40 partially overlapping PCR fragments per individual, which is a time- and resource-consuming process. We have developed a high-throughput, accurate, fast, and cost-effective human MT DNA sequencing approach. In this setup we first generate long-range PCR products for two partially overlapping 7.7 and 9.2 kb MT DNA-specific amplicons, add sample-specific barcodes, and sequence these on the PacBio RSII system to obtain full-length MT DNA sequences for genotyping/haplotyping purposes.

Whole-transcriptome analysis of flow-sorted cervical cancer samples reveals that B cell expressed TCL1A is correlated with improved survival.

Cervical cancer is typically well infiltrated by immune cells. Because of the intricate relationship between cancer cells and immune cells, we aimed to identify both cancer cell and immune cell expressed biomarkers. Using a novel approach, we isolated RNA from flow-sorted viable EpCAM+ tumor epithelial cells and CD45+ tumor-infiltrating immune cells obtained from squamous cell cervical cancer samples (n = 24). Total RNA was sequenced and differential gene expression analysis of the CD45+ immune cell fractions identified TCL1A as a novel marker for predicting improved survival (p = 0.007). This finding was validated using qRT-PCR (p = 0.005) and partially validated using immunohistochemistry (p = 0.083). Importantly, TCL1A was found to be expressed in a subpopulation of B cells (CD3-/CD19+/CD10+/CD34-) using multicolor immunofluorescence. A high TCL1A/CD20 (B cell) ratio, determined in total tumor samples from a separate patient cohort using qRT-PCR (n = 52), was also correlated with improved survival (p = 0.027). This is the first study demonstrating the prognostic value of separating tumor epithelial cells from tumor-infiltrating immune cells and determining their RNA expression profile for identifying putative cancer biomarkers. Our results suggest that intratumoral TCL1A+ B cells are important for controlling cervical cancer development.

MicroRNA classifier and nomogram for metastasis prediction in colon cancer.

BACKGROUND: Colon cancer prognosis and treatment are currently based on a classification system still showing large heterogeneity in clinical outcome, especially in TNM stages II and III. Prognostic biomarkers for metastasis risk are warranted as development of distant recurrent disease mainly accounts for the high lethality rates of colon cancer. miRNAs have been proposed as potential biomarkers for cancer. Furthermore, a verified standard for normalization of the amount of input material in PCR-based relative quantification of miRNA expression is lacking. METHODS: A selection of frozen tumor specimens from two independent patient cohorts with TNM stage II-III microsatellite stable primary adenocarcinomas was used for laser capture microdissection. Next-generation sequencing was performed on small RNAs isolated from colorectal tumors from the Dutch cohort (N = 50). Differential expression analysis, comparing in metastasized and nonmetastasized tumors, identified prognostic miRNAs. Validation was performed on colon tumors from the German cohort (N = 43) using quantitative PCR (qPCR). RESULTS: miR25-3p and miR339-5p were identified and validated as independent prognostic markers and used to construct a multivariate nomogram for metastasis risk prediction. The nomogram showed good probability prediction in validation. In addition, we recommend combination of miR16-5p and miR26a-5p as standard for normalization in qPCR of colon cancer tissue-derived miRNA expression. CONCLUSIONS: In this international study, we identified and validated a miRNA classifier in primary cancers, and propose a nomogram capable of predicting metastasis risk in microsatellite stable TNM stage II-III colon cancer. IMPACT: In conjunction with TNM staging, by means of a nomogram, this miRNA classifier may allow for personalized treatment decisions based on individual tumor characteristics.

Strategies for rapidly mapping proviral integration sites and assessing cardiogenic potential of nascent human induced pluripotent stem cell clones.

Recent methodological advances have improved the ease and efficiency of generating human induced pluripotent stem cells (hiPSCs), but this now typically results in a greater number of hiPSC clones being derived than can be wholly characterized. It is therefore imperative that methods are developed which facilitate rapid selection of hiPSC clones most suited for the downstream research aims. Here we describe a combination of procedures enabling the simultaneous screening of multiple clones to determine their genomic integrity as well as their cardiac differentiation potential within two weeks of the putative reprogrammed colonies initially appearing. By coupling splinkerette-PCR with Ion Torrent sequencing, we could ascertain the number and map the proviral integration sites in lentiviral-reprogrammed hiPSCs. In parallel, we developed an effective cardiac differentiation protocol that generated functional cardiomyocytes within 10 days without requiring line-specific optimization for any of the six independent human pluripotent stem cell lines tested. Finally, to demonstrate the scalable potential of these procedures, we picked 20 nascent iPSC clones and performed these independent assays concurrently. Before the clones required passaging, we were able to identify clones with a single integrated copy of the reprogramming vector and robust cardiac differentiation potential for further analysis.

TSSV: a tool for characterization of complex allelic variants in pure and mixed genomes.

MOTIVATION: Advances in sequencing technologies and computational algorithms have enabled the study of genomic variants to dissect their functional consequence. Despite this unprecedented progress, current tools fail to reliably detect and characterize more complex allelic variants, such as short tandem repeats (STRs). We developed TSSV as an efficient and sensitive tool to specifically profile all allelic variants present in targeted loci. Based on its design, requiring only two short flanking sequences, TSSV can work without the use of a complete reference sequence to reliably profile highly polymorphic, repetitive or uncharacterized regions. RESULTS: We show that TSSV can accurately determine allelic STR structures in mixtures with 10% representation of minor alleles or complex mixtures in which a single STR allele is shared. Furthermore, we show the universal utility of TSSV in two other independent studies: characterizing de novo mutations introduced by transcription activator-like effector nucleases (TALENs) and profiling the noise and systematic errors in an IonTorrent sequencing experiment. TSSV complements the existing tools by aiding the study of highly polymorphic and complex regions and provides a high-resolution map that can be used in a wide range of applications, from personal genomics to forensic analysis and clinical diagnostics. AVAILABILITY AND IMPLEMENTATION: We have implemented TSSV as a Python package that can be installed through the command-line using pip install TSSV command. Its source code and documentation are available at https://pypi.python.org/pypi/tssv and http://www.lgtc.nl/tssv.

Transcriptome and genome sequencing uncovers functional variation in humans.

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project--the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Reproducibility of high-throughput mRNA and small RNA sequencing across laboratories.

RNA sequencing is an increasingly popular technology for genome-wide analysis of transcript sequence and abundance. However, understanding of the sources of technical and interlaboratory variation is still limited. To address this, the GEUVADIS consortium sequenced mRNAs and small RNAs of lymphoblastoid cell lines of 465 individuals in seven sequencing centers, with a large number of replicates. The variation between laboratories appeared to be considerably smaller than the already limited biological variation. Laboratory effects were mainly seen in differences in insert size and GC content and could be adequately corrected for. In small-RNA sequencing, the microRNA (miRNA) content differed widely between samples owing to competitive sequencing of rRNA fragments. This did not affect relative quantification of miRNAs. We conclude that distributing RNA sequencing among different laboratories is feasible, given proper standardization and randomization procedures. We provide a set of quality measures and guidelines for assessing technical biases in RNA-seq data.

Deep sequencing of RNA from immune cell-derived vesicles uncovers the selective incorporation of small non-coding RNA biotypes with potential regulatory functions.

Cells release RNA-carrying vesicles and membrane-free RNA/protein complexes into the extracellular milieu. Horizontal vesicle-mediated transfer of such shuttle RNA between cells allows dissemination of genetically encoded messages, which may modify the function of target cells. Other studies used array analysis to establish the presence of microRNAs and mRNA in cell-derived vesicles from many sources. Here, we used an unbiased approach by deep sequencing of small RNA released by immune cells. We found a large variety of small non-coding RNA species representing pervasive transcripts or RNA cleavage products overlapping with protein coding regions, repeat sequences or structural RNAs. Many of these RNAs were enriched relative to cellular RNA, indicating that cells destine specific RNAs for extracellular release. Among the most abundant small RNAs in shuttle RNA were sequences derived from vault RNA, Y-RNA and specific tRNAs. Many of the highly abundant small non-coding transcripts in shuttle RNA are evolutionary well-conserved and have previously been associated to gene regulatory functions. These findings allude to a wider range of biological effects that could be mediated by shuttle RNA than previously expected. Moreover, the data present leads for unraveling how cells modify the function of other cells via transfer of specific non-coding RNA species.

New methods for next generation sequencing based microRNA expression profiling.

BACKGROUND: MicroRNAs are small non-coding RNA transcripts that regulate post-transcriptional gene expression. The millions of short sequence reads generated by next generation sequencing technologies make this technique explicitly suitable for profiling of known and novel microRNAs. A modification to the small-RNA expression kit (SREK, Ambion) library preparation method for the SOLiD sequencing platform is described to generate microRNA sequencing libraries that are compatible with the Illumina Genome Analyzer. RESULTS: High quality sequencing libraries can successfully be prepared from as little as 100 ng small RNA enriched RNA. An easy to use perl-based analysis pipeline called E-miR was developed to handle the sequencing data in several automated steps including data format conversion, 3' adapter removal, genome alignment and annotation to non-coding RNA transcripts. The sample preparation and E-miR pipeline were used to identify 37 cardiac enriched microRNAs in stage 16 chicken embryos. Isomir expression profiles between the heart and embryo were highly correlated for all miRNAs suggesting that tissue or cell specific miRNA modifications do not occur. CONCLUSIONS: In conclusion, our alternative sample preparation method can successfully be applied to generate high quality miRNA sequencing libraries for the Illumina genome analyzer.

Comprehensive gene-expression survey identifies wif1 as a modulator of cardiomyocyte differentiation.

During chicken cardiac development the proepicardium (PE) forms the epicardium (Epi), which contributes to several non-myocardial lineages within the heart. In contrast to Epi-explant cultures, PE explants can differentiate into a cardiomyocyte phenotype. By temporal microarray expression profiles of PE-explant cultures and maturing Epi cells, we identified genes specifically associated with differentiation towards either of these lineages and genes that are associated with the Epi-lineage restriction. We found a central role for Wnt signaling in the determination of the different cell lineages. Immunofluorescent staining after recombinant-protein incubation in PE-explant cultures indicated that the early upregulated Wnt inhibitory factor-1 (Wif1), stimulates cardiomyocyte differentiation in a similar manner as Wnt stimulation. Concordingly, in the mouse pluripotent embryogenic carcinoma cell line p19cl6, early and late Wif1 exposure enhances and attenuates differentiation, respectively. In ovo exposure of the HH12 chicken embryonic heart to Wif1 increases the Tbx18-positive cardiac progenitor pool. These data indicate that Wif1 enhances cardiomyogenesis.